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ABSTRACT 


A  theoretical  solution  is  presented  for  species  diffusion  in  the  com¬ 
pressible  laminar  boundary  layer  on  a  flat  plate  for  frozen  gas-phase 
chemistry.  It  is  assumed  that  the  recombination  at  the  wall  is  a  first-order 
reaction,  and  the  chordwise  concentration  gradients  are  neglected.  The 
solution  parallels  Crocco's  treatment  except  for  the  assumption  governing 
the  viscosity  variations  with  temperature.  A  modified  Chapman-Rube  sin 
constant  is  used  here,  and  it  is  assumed  that  the  product  of  the  density  and 
viscosity  varies  as  a  power  of  the  static  enthalpy-  This  makes  it  possible 
to  closely  approximate  the  actual  yO/t- product  in  the  vicinity  of  the  wall. 

The  solution  for  the  species  equation  and  the  energy  equation  are 
obtained  in  terms  of  the  shear  in  the  boundary  layer.  The  latter  is  deter¬ 
mined  following  an  empirical  calculation  used  by  Young.  The  result  given 
here  should  be  somewhat  more  general  since  the  Frandtl  number  effects 
are  included  in  a  more  complete  fashion.  The  solutions  for  heat  transfer 
and  skin  friction  are  compared  with  existing  exact  numerical  solutions  to 
show  that  they  agree  with  typical  errors  of  about  5  percent. 
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I.  INTRODUCTION 


A  program  of  research  is  in  progress  at  CAL  aimed  at  developing  a 
catalytic  probe  as  a  diagnostic  device  for  measuring  species  concentrations 
in  nozzle  expansions  of  high  temperature  air.  This  device  was  originally 
suggested  by  Hoenig,^  and  has  received  increasing  attention  in  recent  years. 
Basically,  it  consists  of  two  adjacent  heat  transfer  gages.  One  is  coated 
with  a  surface  which  is  catalytic  to  atom  recombination  while  the  other  is 
coated  with  a  noncatalytic  surface.  The  difference  between  the  heat  transfer 
rates  sensed  by  the  two  gages  is  proportional  to  the  number  of  atoms  striking 
the  surface,  and  this  can  be  related  to  the  ambient  atom  concentration  if  an 
appropriate  theoretical  solution  is  available. 

Theoretical  studies  of  catalytic  surfaces  in  reacting  flows  have  been 

pursued  for  about  a  decade  and  solutions  are  presently  available  for  species 

2-11 

diffusion  in  the  laminar  boundary  layer  both  for  the  frozen  case  and  for 

12-14 

the  case  where  there  is  species  production.  A  variety  of  configurations 

are  treated  in  these  references  including  the  flat  plate  with  and  without 
pressure  gradient,  wedges,  cones,  and  blunt  bodies.  These  also  cover  a 
range  of  assumptions  for  the  transport  properties  including  a  linear  viscosity- 
temperature  relation,  constant  Prandtl  and  Schmidt  numbers  which  differ 
from  unity,  and  variable  Prandtl  and  Schmidt  numbers.  With  this  wide 
assortment  of  theoretical  solutions  available,  it  was  expected  that  experi¬ 
mental  data  could  be  interpreted  within  the  framework  of  at  least  one  of  them. 
Early  experimental  results  showed,  however,  that  there  were  important 
effects  apparently  stemming  from  the  transport  parameters  {Prandtl  number 
and  Schmidt  number)  as  well  as  variations  in  the  product  of  the  density  and 

viscosity.  This  might  have  been  anticipated  from  the  numerical  solutions 

g 

given  by  Janowitz  and  Libby.  In  that  work,  numerical  solutions  were  obtained 
for  the  frozen  laminar  boundary  layer  for  certain  specific  cases.  Those 
included  constant  transport  parameters  as  well  as  realistic  models  for  the 
transport  processes  in  the  boundary  layer.  The  comparison  of  the  computed 
results  obtained  with  the  usual  assumptions  for  the  transport  parameters 
{cr  =  J>c=  1 ,  =  constant)  discloses  large  and  important  effects  stemming 
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from  these  approximations.  The  effects  of  Prandtl  and  Schmidt  numbers 

5 

that  differ  from  unity  can  be  accounted  for  with  Inger's  extension  of  Lees' 

1  5 

theory.  However,  all  of  the  cited  solutions  use  the  assumption, 

S>M- 

—  =  constant.  The  purpose  of  the  present  report  is  to  present  an 

'c 

approximate  theoretical  solution  for  species  diffusion  in  the  compressible 
laminar  boundary  layer  on  a  flat  plate  which  accounts  for  Prandtl  and 
Schmidt  numbers  that  differ  from  unity  and  allows  for  a  p/L-  product  that 
varies  as  a  power  of  the  static  enthalpy. 

A  powerful  approach  in  treating  the  species  diffusion  in  the  laminar 

boundary  layer  is  to  take  advantage  of  the  fact  that  with  the  boundary  layer 

assumptions,  the  species  conservation  equation  is  decoupled  from  the 

momentum  and  energy  equations  if  it  is  assumed  that  the  Lewis  number  is 

unity.  This  makes  it  possible  to  extend  existing  solutions  by  writing  down 

the  species  conservation  equation  and  solving  it  in  terms  of  the  solution  for 

either  the  momentum  or  the  energy  equation.  That  approach  is  used  here 

and  is  applied  to  Crocco's  solution  for  the  flat  plate. Crocco's  solution, 
17  18 

as  reported  and  extended  by  Young,  is  used  here  because  it  covers  the 
cases  of  interest,  and  because  the  results  for  skin  friction  and  heat  transfer 
are  expressible  in  simple  yet  accurate  semiempirical  formulas. 

In  the  succeeding  paragraphs,  Crocco's  solutionis  extended  to  include 
species  diffusion  and  a  solution  for  the  species  equation  in  terms  of  the 
solution  to  the  momentum  equation.  Following  this,  Young's  semiempirical 
generalization  of  Crocco's  results  for  the  skin  friction  and  heat  transfer  is 
reviewed  and  extended  to  better  account  for  Prandtl  number  effects.  The 
model  for  the  p jjl -variations  in  the  boundary  layer  is  critically  examined 
next  to  arrive  at  meaningful  method  for  establishing  the  free  constants. 
Finally,  the  present  results  are  compared  with  exact  numerical  solutions 
to  assess  their  accuracy. 
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II.  THEORETICAL  DEVELOPMENT 


1  G  17 

The  succeeding  derivation  follows  that  of  Crocco,  *  and  his 
development  is  repeated  here  in  the  interests  of  completeness.  Since  the 
original  theory  was  for  an  ideal  gas,  first  it  must  be  shown  that  the  energy 
equation  for  a  nonreacting  dissociated  gas  can  be  reduced  to  the  energy 

19 

equation  used  by  Crocco.  Using  the  energy  equation  given  by  Dorrance, 
the  procedure  consists  of  subtracting  out  the  momentum  equation  and 
defining  a  mean  specific  heat  for  the  gas  mixture, 


cp  -  C 

where  oc.  and  C  are  the  mass  fraction  and  the  specific  heat  at  constant 
i  A  c 

pressure,  respectively,  for  the  i  species.  The  static  enthalpy  is  then 


/>  =  E  .  h  -  I  -r  L 7  «  -  /> 

L  L  t  L  Di 

where  I  -  CpT  and  hD_  is  the  heat  of  formation  for  the  i  ih  species.  Using 
this  in  the  above  form  of  the  energy  equation,  the  final  form  of  the  energy 
equation  is  obtained  after  subtracting  out  the  species  conservation  equation. 
Assuming  zero  pressure  gradient,  a  constant  Prandtl  number,  cr'  ,  and  a 
unit  Lewis  number,  the  system  of  governing  relations  is 


momentum 

continuity 
energy 

species 


dn.  dn 

PA  nrr-  ^P^ 


3X 


dy 


A 


dU. 

d<j 


dx 


(p'A-)+Jj[S>v')~  0 


di  ,  dl 

P^—  +  P^  — —  =  /* 


3tC 


(  a  u 

(  di 


3  oc.  3°c  - 

/OV, - -  +  JO  v  - —  = 

y  dX  dij 


1 


Sc 


cr  dtj 
d<j  (  dy 


,  di  ) 
^  3ij  ) 


with  the  boundary  conditions 


for 

n 

O 

U.-0, 

n 

b* 

i 

1! 

for 

ij  =  <z>  ; 

u— , 
c’ 

I-*l 

& 

A* 


oc  ■ 


(1) 

(2) 

(3) 

{4} 


P*  (Sc  Dy 
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Equations  1  through  3  and  the  accompanying  boundary  conditions  are 
identically  those  treated  by  Crocco.  Equation  4  is  the  additional  equation 
to  be  solved  subject  to  the  indicated  first-order  wall  boundary  condition. 


Following  the  Crocco  method  of  solution,  the  independent  variables 


are  changed  from  {  x  ,  u  )  to  (  x  ,  ic  )  and  a  new  dependent  variable,?^  a  , 
"  3  y 

is  introduced.  Making  this  transformation  and  then  combining  the  resulting 
momentum  and  continuity  equations,  the  transformed  system  of  equations  is 


momentum  and  continuity 

u — 

«  0 

(5a) 

r  -|  3T  dl 

?  u 

z 

- 

dx  L  zr  J 

3'1 1 

dr 

(6a) 

energy  Q  -  <r]  r  —  — 

+  r 

<r  t 

- -  -  CT/O  JX  K  -  =  0 

du.  d  x 

dx<x.  r  i 

species  Z  ^  a,  +  ^ 

„  dz 
z  — 

d  ix 

doi^ 

3lc 

-  Sc/OyU  u. 

doc.  L 

L  =  0 

dX 

(7a) 

with  the  boundary  conditions 

v-  =•  o  I  =  iw ,  zr  =  rw ,  (o,  x.)  = 


5, 


—  U^-  L  —  Xf  ,  Z  —  O ,  oc  ^ 

Crocco's  analysis  anticipates  similar  velocity  profiles  and  enthalpy 
profiles,  and  therefore,  neglects  the  chordwise  variations  in  these  quantities. 
This  step  is  consistent  with  the  zero  pressure  gradient  assumption  if  it  is 
assumed  that  the  wall  temperature  is  invariant.  It  will  be  assumed  here 
that  the  chordwise  variations  in  species  concentrations  are  also  negligible. 
This  is  consistent  with  the  other  assumptions  and  is  valid  if  the  wall  con¬ 
centrations  are  small  compared  with  those  in  the  ambient  stream  (the 
diffusion-controlled  regime),  or  if  the  wall  concentration  approaches  the 
ambient  concentration  (the  rate-controlled  regime).  These  regimes  will  be 
defined  in  a  more  quantitative  manner  after  obtaining  the  solution.  Antici¬ 
pating  similar  profiles,  it  is  assumed  that  the  shear  is  of  the  form 

rr  =  0  (x)  f  (u,) 
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which  when  substituted  in  Equation  5a  yield  4>  = 
become 


1 


Equation  5a  through  7a 


dzf 

d  it x 


duz 


(5b) 

(6b) 

{7b} 


Prior  to  solving  the  above  relations,  they  are  made  dimensionless  by  dividing 
by  the  quantities  at  the  edge  of  the  boundary  layer 


w==*<r? 


7=2-^ 


ot.  -  Ot;  L 
i  Le 


and  by  introducing  the  new  variable 

*£(?)  =  fM  j 


S>e^e  ue 


(8) 


(9) 


Until  now,  we  have  duplicated  Crocco's  analysis  in  every  detail.  At 
this  point  it  is  necessary  to  depart  from  that  procedure  in  one  respect.  If 
the  relations  in  Equation  8  are  substituted  in  Equation  5b,  the  second  term 
becomes  j  ??  ,  and  Crocco  then  makes  the  assumption  so  that 


for  an  ideal  gas,  the  momentum  equation  becomes 


It  is  not  possible  to  proceed  as  directly  in  the  present  problem  because  the 
species  concentrations  vary  in  the  boundary  iayer  and  consequently,  influence 
the  viscosity  and  the  molecular  weight.  Instead,  it  is  observed  that  the 
power-law  assumption  for  the  temperature-viscosity  relation  is  made 
primarily  for  the  purposes  of  reducing  the  y?/<.-product  to  a  simple  function 
of  the  enthalpy.  Consequently,  we  will  proceed  directly  to  an  equivalent 
assumption. 


/A 
/°e  /*t 


(10) 
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and  will  defer  a  critical  appraisal  of  this  assumption  until  a  later  section. 

It  should  be  noted,  however,  that  C'  is  a  constant  analogous  to  the 

20  * 

Chapman-Rube  sin  constant  or  to  Cheng's  improvement  on  that  constant. 
Using  the  dimensionless  quantities  in  Equation  8,  the  assumption  given  by 
Equation  10,  and  the  dimensionless  variable, 


^1?)  =  /  (u) 


(H) 


the  final  equations  are  as  follows 


F  7T*  * 


7 


**  w ** 


rff 


-f 


=  0 


Id? 


ut 

2  +  ^  I 


e  j 


di  o(F 


F - \  + 


r  1  dF  d  £ 

L'-^J  77  Fz=° 


~  o 


(5c) 

(6c) 

(7c) 


with  the  boundary  conditions 

C./& 

7  =  0 :  F  «  F(o)  =  —=7  ,  4=  Ko) 


7^, 

C-f>e  7^ 


7  =  O  :  F  =  O,  t'  =  T=  1 

Equations  5c  and  6c  are  identical  to  Crocco's  governing  equations  and  his 
solutions  apply  directly  if  it  is  noted  that  the  nondimensional  shear  parameter. 
Equation  11,  differs  from  that  used  by  Crocco,  Crocco  numerically  integrated 
Equation  5c  to  obtain  solutions  to  the  momentum  equation.  Both  the  energy 
equation  and  the  species  equation  can  be  integrated  in  terms  of  the  non- 
dimensional  shear.  The  solution  for  the  enthalpy  distribution  is  that  given 
in  Reference  17 


c  (7)  =  i  (O)  + 


A,  (7,  O')  -  (T 


(12) 


whe  re 


Ff%) 

F(O) 


1-<r 


drl„  i  d7l, 
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are  tabulated  in  Reference  17  as  a  function  of  ??  and  cr'  .  The  values  of 
these  quantities  evaluated  at  the  edge  of  the  boundary  layer  enter  into  the 
solution  and  are  accurately  represented  (errors  of  less  than  2  percent)  by 

At  (  f,  tr)  =  o'  ~ 1/3  2  Bt  (>,  <r)  =  ,/x- 


The  species  conservation  equation,  Equation  7c,  can  be  integrated  immedi- 

i 

ately,  and  after  applying  the  boundary  conditions,  the  solution  is 


)  = 


L 

£  (zio') 

1  i-E(Z,(r) 


■h 


1  +  £  (x,<r) 


(13) 


where 


EM  = 


Cf  /~Rt 

A ,  (/, 

J  Z<t/r& 


'  Cf/*2 

1 

P  K 

_  2  dr 

The  similarity  assumption,  neglecting  the  chordwise  variations  in 
species  concentrations,  can  now  be  assessed.  Since  the  present  analysis 
is  restricted  to  zero  pressure  gradient,  the  only  chordwise  variations  must 
result  from  variations  in  the  wall  concentrations.  Equation  13  shows  that 
the  wall  concentration  is 


£  (x,r) 
f  +  £  (x,  (T) 


/nr  £  (z,#")  — ►  0 
Ar  £  (*,o -’)—+■  a> 


Consequently,  the  similarity  assumption  is  valid  in  the  diffusion-controlled 
regime  where  £(%,<r)  is  vanishingly  small,  or  equivalently,  at  large  Reynolds 
numbers.  The  assumption  is  also  valid  in  the  rate-controlled  regime,  where 
the  wall  reaction  speed,  ,  is  vanishingly  small,  or  equivalently,  where 
the  catalytic  efficiency,  y  ,  approaches  zero. 
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III.  THE  SURFACE  HEAT  TRANSFER  AND  SKIN  FRICTION 


The  heat  transfer  to  a  surface  in  a  gas  with  a  single  reacting  species 


,  „  19 

is  given  by  Dorrance  as 


dl 


Lc, 


- p 

Defining  a  Stanton  number  as  C  = 


oc 


d  (j  J 


(14) 


-s 


«e(“r-  r„) 


■where  H  =/  +  -  uzt  «.  h. 


and  using  the  solutions  as  given  by  Equations  12  and  13,  the  heat  transfer  to 
the  wall  (  1  )  can  be  expressed  in  terms  of  the  skin  friction  coefficient 

as  follows. 


Cr  /fit 


z(or) 


2/s 


1  r 


</z 


[H 


Od  k 
£  nQ 


w 


(15) 


where 

4  =  *’ 

"  1  *  E(X,<r) 

The  skin  friction  was  determined  by  Crocco  by  numerically  integrating 

1 8 

the  momentum  equation.  Young  subsequently  generalized  those  results  in 
a  semiempirical  calculation  to  obtain  a  closed-form  expression  for  the  skin 
friction.  Young's  calculation  will  be  extended  here  to  better  account  for  the 
influence  of  the  Prandtl  number.  The  basis  for  the  calculation  follows  from 
Equation  5c  which  can  be  put  in  the  following  form 


F(o)  = 


£1  IQ 


(16) 


Young  notes  that  Crocco1  s  numerical  results  for  F(y)  are  closely  approximated 
by  F(0)J f-T?3-  which,  substituted  in  Equation  16  yields  an  approximate  expression 
for  the  skin  friction  coefficient. 


(17) 
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Equation  17  can  be  integrated  if  the  enthalpy  distribution  in  the  boundary 
layer  is  known.  Young  uses  the  enthalpy  for  a  unit  Prandtl  number  and 
introduces  a  Prandtl  number  dependence  after  completing  the  integration 
by  adjusting  the  final  expression  to  agree  with  known  exact  solutions.  The 
approach  used  here  is  to  retain  Young's  expression  for  the  shear  distribution 
in  the  boundary  layer,  but  to  fit  a  power  series  to  Crocco's  solution  for  the 
enthalpy  distribution,  or  equivalently,  a  power  series  for  Ar (??■>&')  and  8f{ y ><?')• 
The  expressions  used  are 

Ai(rl'<r)  —  7  a,  ( t r +  °-z  (&*)  Y5 

( 1 8) 

B,  (7?’°')~  ~j  [yz+h,  ] 

The  parameters,  a  >  a  ,  b  ,  and  b  are  plotted  as  a  function  of  the  Prandtl 

T  2  J  £ 

number  in  Figure  1.  The  worst  errors  between  the  results  tabulated  in 
Reference  17  and  Equation  18  were  for  a  Prandtl  number  of  1/2.  Those 
errors  were  about  2  to  3  percent  and  occurred  near  the  edge  of  the  boundary 
layer. 

Equation  17  has  been  integrated  using  Equations  12  and  18  for  the 
enthalpy  distribution.  The  procedure  for  this  integration  was  identical  to 
that  used  by  Young;  namely,  to  note  that  the  expression  for  the  enthalpy  is 
of  the  following  form 


and  that  each  of  the  quantities  in  the  brace  are  always  less  than  unity.  This 
suggests  that  the  enthalpy  function  can  be  approximated  by  the  first  term  in 
the  series  expansion, 
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After  the  integration  has  been  completed,  the  procedure  ia  reversed  to  yield 
the  original  power  dependence.  The  details  of  this  procedure  can  be  obtained 
in  Reference  18.  The  integration  and  subsequent  return  to  the  original  form 
yields  the  following  result 


F(o)  =  Cf  [>-  -Z-]  '  |  Cf  +  CA  +  — —  O.trSl  +  0.04*4  a,(r)-0. 2094 


c*  i 

* 

1_  -f  J  u 


0777  b  0.7379  h  (<r)  -0.0703  b  (<r) 


“I  X 

'J 


It  is  obvious  that  this  procedure  is  approximate  and  can  only  be  justified  by 

1 8 

comparing  the  result  with  exact  solutions.  As  shown  by  Young,  it  does 
however  yield  a  result  in  close  agreement  with  the  Blasius  limit.  For  the 
case  where  ^^7,  F (o)~  =■  0,055  (since  =  1  )  as  compared  with 

0.  664  for  the  Blasius  limit.  This  close  agreement  suggests  the  final  form 
for  the  skin  friction  coefficient. 


Cf/fi2  =  0665  U  (cr)  +  [ /. -  k{  (erf\  i  (o)  +  K  J+) 


(19a) 


The  constants,  kf  (<r)  and  (cr) ,  are  related  to  the  constants  a^,  6  ,  and 
bz  and  are  functions  only  of  the  Prandtl  number.  They  are  plotted  in 
Figure  2  and  compared  with  the  values  recommended  by  Young.  It  can  be 
seen  that  they  are  in  reasonable  agreement  at  cr"-  1  ,  but  the  parameters 
arrived  at  here  show  a  rather  strong  Prandtl  number  dependence.  They  can 
be  approximated  by  the  following  relations 

Kt  (<r)  ^  0.55  er^7  (cr)  ^  0-0965  jjr" 


-Vs 


+  cr 


'°/3 


to  yield  for  the  skin  friction 


*  t&Il 

Cj,  fRc  ~  0.665  fc^i^L  (0)  -b  0.55<rh  \j~  i(0f\  i-  o.  o  9  65  ]£]  *  (19b) 


Equations  15  and  19  are  the  final  relations  for  the  heat  transfer  and 

skin  friction  coefficients.  They  both  depend  on  two  constants,  c'  and  cO 

*  eff 

which  are  as  yet  unspecified.  These  will  be  determined  in  the  next  section 
where  consideration  is  given  to  the  significance  of  these  parameters. 
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IV.  DENSITY-VISCOSITY  PRODUCT 


A  basic  difficulty  in  any  compressible  laminar  boundary  layer  develop¬ 
ment  is  in  treating  the  product  of  the  density  and  the  viscosity.  There  is 
a  variety  of  approximate  techniques  used  in  the  theory  for  nonreacting  gases, 
and  it  is  natural  to  look  to  that  literature  for  a  guide.  The  simplest  approach 
in  nonreacting  gases  is  to  assume  that  the  viscosity  is  linearly  related  to  the 
temperature  so  that  within  the  boundary  layer  approximations,  the  />/>■- 
product  is  constant  through  the  boundary  layer.  This  procedure  results  in 
the  familiar  Chapman-Rubesin  constant,  which  makes  it  possible  to  match 

the  viscosity  at  one  point  in  the  flow  field.  The  original  choice^  was  at  the 

2 1 

wall,  and  it  was  found  to  be  accurate  at  supersonic  speeds.  Cheng  refined 

this  procedure  by  matching  the  viscosity  at  the  point  in  the  boundary  layer 

corresponding  to  the  average  temperature.  This  approach  has  been  used  at 

relatively  high  temperatures  (3000°K)  and  shown  to  produce  accurate  results 

2122 

when  compared  with  experimental  data.  >  It  has  been  shown  to  lead  to 

23 

serious  errors,  however,  when  applied  at  temperatures  of  5000  to  6000  “K. 
The  power-law  approximation  used  by  Crocco  is  useful  since  it  makes  it 
possible  to  match  the  viscosity  at  one  point  in  the  flow  field  and  to  approxi¬ 
mate  the  actual  temperature  dependence.  A  better  approximation  is  to 
introduce  the  Chapman-Rubesin  constant  in  conjunction  with  the  power-law 
approximation  as  done  in  Reference  23.  This  makes  it  possible  to  match  the 
viscosity  at  two  points  in  the  flow  field  while  approximating  the  actual 
temperature  dependence. 

None  of  these  approximations  to  the  viscosity  is  useful  in  the  present 
problem  owing  to  the  complications  introduced  by  the  variable  species  con¬ 
centrations  in  the  boundary  layer.  The  species  concentration  changes  the 
molecular  weight  and  thereby  introduces  a  dependence  on  the  concentrations 
as  well  as  on  the  temperature..  Depending  on  the  model  used  for  estimating 
the  viscosity,  the  species  concentrations  may  or  may  not  enter  into  that 
quantity  too.  Because  of  these  complications,  the  approach  used  here  is  to 
relate  the  /5/t-product  to  the  enthalpy  using  an  approximation  analogous  to 
that  used  in  approximating  just  the  viscosity. 
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The  approximation  introduced  in  Equation  10  is  analogous  to  that  used 
in  Reference  23 


-  W 


cj  - t 
Hf 


except  that  it  must  be  noted  that  the  exponent,  ,  is  not  necessarily  the 

exponent  normally  used  in  relating  the  viscosity  to  the  temperature.  It  and 
C*  are  free  constants  to  be  determined  by  matching  the  actual  /)/*■  -product 
at  two  points  in  the  boundary  layer.  The  two  points  chosen  here  are  at  the 
wall  and  at  the  point  in  the  boundary  layer  corresponding  to  the  average  value 
of  C.  This  procedure  is  analogous1  to  that  used  in  References  21  and  23, 
where  the  average  temperature  is  used  to  determine  the  modified  Chapman- 
Rubesin  constant.  The  second  point  is  chosen  at  the  wall  because  interest 
centers  on  the  transport  processes  at  the  wall. 

The  average  value  of  the  enthalpy,  i  ,  can  be  determined  by  integrating 
Equation  12,  using  Equation  18  to  approximate  the  Prandtl  number  dependence. 
The  result  is 

3. 

i  =4  (o)  +  *r3(r)  t  k*  «r>  jjS-  (20a) 

The  constants  k3(o')  and  K^(cr)  are  functions  only  of  the  Prandtl  number  and 
are  plotted  in  Figure  3.  They  can  be  closely  approximated  by  the  relations, 

K  (ff~)  =5=  ~  <rZ/l  K  (<r-)  m  ~  o'  /S 

5  2  4  6 


and  the  relation  for  the  average  enthalpy  becomes 

£  =  £(0)  +  -~o^7  [f-  if oj]  +  j  er *'5  (20b) 

*  6 

The  average  enthalpy  and  the  species  concentration  corresponding  to  the 
average  value  of  the  enthalpy,  £  U  ),  are  used  to  determine  the /> ju.- product 
that  corresponds  to  the  average  value  of  the  enthalpy.  The  species  concen¬ 
tration  is  computed  from  Equation  13  using  the  value  of  corresponding  to  the 
average  enthalpy,  r^( i )  .  For  unit  Prandtl  number,  this  value  can  be  deter¬ 
mined  from  the  quadratic  relation, 


re 


h  *  (<r)  *  k3  (er) 


1-  i  (0) 


=  0 


(21) 
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For  Prandtl  numbers  other  than  unity,  the  value  of  can  be  determined 
from  Figure  4.  The  free  constants  and  are  then  given  by  the 

following  relations. 


J07? 

1  -  - 

/, T7 


ygv  >v 
_/)(!)  /4(t)_ 

i 

_  i  (0)  _ 


(22) 


S>e  /*e 
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V.  COMPARISON  WITH  EXACT  SOLUTIONS 


The  present  theory  has  been  applied  to  those  cases  treated  by 

g 

Janowitz  and  Libby,  and  the  results  are  compared  with  the  exact  numerical 

solutions  in  Table  I.  The  present  theory  was  evaluated  using  the  gas  model 

and  the  relations  for  the  transport  properties  that  were  used  in  Reference  8. 

The  only  additional  computation  made  here  was  to  calculate  the  total 

temperature  at  the  edge  of  the  boundary  layer  using  real  gas  charts  for 

nitrogen  and  taking  the  density  altitude  as  250,  000  feet.  The  theory  of 

Reference  8  makes  the  same  similarity  assumptions  used  here.  Consequently, 

* 

the  comparison  sheds  no  light  on  the  accuracy  of  that  approximation. 

There  are  four  basic  flows  covered  in  Table  I.  The  first  is  one  in 
which  the  Prandtl  and  Schmidt  numbers  are  unity  and  the  />M~  ratio  is  equal 
to  one.  This  is  the  Blasius  limit,  and  Table  I  shows  that  the  present 
treatment  reduces  to  that  limit  with  good  accuracy. 

The  second  basic  flow  covered  in  Table  I  is  one  in  which  the  /?/<-ratio 
is  unity,  and  the  Prandtl  and  Schmidt  numbers  are  0.7,  These  cases  provide 
a  check  on  the  methods  used  to  account  for  the  transport  parameters  in 
computing  the  heat  transfer,  and  Table  I  shows  that  the  solutions  agree 
within  3  percent. 

The  third  flow  situation  in  Table  I  is  one  in  which  the  transport 
parameters  are  unity,  but  the  /} jjl -ratio  varies  through  the  boundary  layer. 
This  flow  provides  an  excellent  test  of  the  present  theory  since  it  closely 
fits  the  theoretical  assumptions  {&.  =  1,  5C  =  o' =  constant).  Consequently, 
the  comparison  will  furnish  a  good  check  on  the  method  used  to  account  for 
the  /3/i- variation  in  the  boundary  layer.  The  calculations  for  this  case 
followed  the  procedures  outlined  in  the  previous  sections  with  one  addition. 


The  accuracy  of  the  similarity  assumption  can  be  assessed  from  the  results 
given  by  Chambre  and  Acrivos.  They  test  it  over  the  complete  range, 
from  zero  to  infinite  catalycity,  and  show  that  the  largest  error  in  wall 
concentration  is  about  10  percent  occurring  at  values  of  E(x,r)ot  about  1/2. 
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In  determining  the  value  of  /?/<  corresponding  to  the  average  enthalpy,  i  , 
the  quadratic  relation  for  Equation  21,  yielded  two  values.  This  is 

due  to  the  fact  that  the  wall  temperature  and  the  temperature  at  the  edge  of 
the  boundary  layer  are  nearly  equal.  The  value  used  was  that  closest  to 
the  wall.  This  procedure  is  consistent  with  the  objective  of  matching  the 
/>/* -ratio  in  the  vicinity  of  the  wall.  The  coinparison  -with  the  exact 
solutions  in  Table  I  shows  that  the  present  theory  predicts  the  skin  friction 
with  errors  of  less  than  5  percent.  The  two  solutions  for  heat  transfer 
show  typical  discrepancies  of  3  to  4  percent,  though  in  one  instance  (Case  2) 
the  discrepancy  is  7  percent.  In  general,  it  appears  that  the  theory  is  quite 
accurate  when  applied  to  these  cases. 


The  fourth  situation  covered  in  Table  I  closely  corresponds  to  an 
actual  flow  in  that  the  transport  parameters  and  the  />/*■- ratio  vary  through 
the  boundary  layer.  The  calculations  for  this  case  should  follow  an  iteration 
procedure  to  determine  the  atom  concentrations  and  the  transport  parameters 
at  the  wall.  Instead,  the  transport  parameters  and  the  concentrations  at  the 
wall  given  by  Janowitz  and  Libby  were  used  to  calculate  the  flow  properties 
used  in  determining  C+  and  *7,^  (Equation  22).  With  these  so  determined, 
the  skin  friction,  heat  transfer,  atom  concentration,  and  /) -ratio  at  the 
wall  were  calculated.  Table  I  shows  that  the  discrepancies  between  the  two 
solutions  for  skin  friction  are  typically  5  percent  and  are  always  less  than 
10  percent.  The  difference  between  the  two  solutions  for  heat  transfer  are 
typically  5  percent,  though  in  one  instance  (Case  8)  the  discrepancy  is 
14  percent  and  in  one  other  (Case  2)  it  is  10  percent.  These  two  errors 
appear  to  stem  largely  from  the  assumption  that  the  Lewis  number  is  unity 
and  can  be  partially  accounted  for  by  using  the  complete  expression  for  heat 
transfer. 


cj_  { 


/z 

Ht~ 


—  R 

H-T  L? 


/- Jfi. 
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If  the  above  expression  is  used  to  compute  the  heat  transfer  for  these  two 
cases,  the  discrepancies  are  less  than  5  percent.  There  is,  of  course,  no 
firm  justification  for  this  procedure  and  it  can  only  be  viewed  as  indicative 
of  the  source  of  the  discrepancies. 


15 


AEOC-TR-67-88 


VI.  CONCLUDING  REMARKS 


Crocco’s  theory  for  the  compressible  laminar  boundary  layer  on  a 
flat  plate  has  been  extended  to  obtain  an  approximate  analytic  boundary  layer 
solution  for  species  diffusion  to  a  catalytic  surface.  This  extension  has 
paralleled  the  development  as  reported  by  Young  except  for  the  viscosity 
assumption  and  the  empirical  treatment  for  the  skin  friction.  Instead  of 
using  Crocco's  assumption  that  the  viscosity  varies  as  a  power  of  the 
temperature,  the  assumption  made  here  was  that  the  product  of  the  density 
and  viscosity  varied  as  a  power  of  the  static  enthalpy  {  C  T).  In  addition, 
a  modified  Chapman-Rubesin  constant  is  introduced  to  improve  the  accuracy 
of  the  assumption.  This  procedure  is  in  the  spirit  of  the  usual  viscosity 
assumption  except  that  it  lumps  together  the  viscosity  variations  and  species 
variations,  as  they  influence  the  -ratio. 

The  empirical  treatment  of  the  skin  friction  given  here  differs  from 
that  jgiven  by  Young  in  that  his  is  based  on  the  enthalpy  profile  for  unit 
Prandtl  number.  The  treatment  given  here  attempts  to  better  account  for 
the  Prandtl  number  variations  by  fitting  a  power  series  to  the  enthalpy 
profiles. 

Solutions  obtained  with  the  present  theory  have  been  compared  with 
exact  numerical  solutions  to  show  that  typically  they  agree  within  about 
5  percent.  This  agreement  is  regarded  as  satisfactory  since  the  theory 
accounts  for  variations  in  the  y?/<. -product  and  since  it  provides  a  relatively 
rapid  method  for  making  estimates  of  the  skin  friction  and  heat  transfer. 
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Figure  3  Solution  for  the  Averoge  Enthlopy 
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